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Abstract 

We present a model treating the kinetics of adsorption of soluble surface-active 
molecules at the interface between an aqueous solution and another fluid phase. The 
model accounts for both the diffusive transport inside the solution and the kinetics 
taking place at the interface using a free-energy formulation. In addition, it offers 
a general method of calculating dynamic surface tensions. Non-ionic surfactants 
are shown, in general, to undergo a diffusion-limited adsorption, in accord with 
experimental findings. 
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Surface-active agents (surfactants) play a major role in various fields, including petro- 
chemical technologies, detergents, biological systems, etc. In some important cases, equi- 
librium properties of the surfactant adsorption are not sufficient, and knowledge of kinetic 
processes is required. In particular, we mention processes of fast wetting, foaming and 
the stability of liquid films. The kinetics of surfactant adsorption have been addressed 
by numerous experimental and theoretical studies, and various experimental techniques 
have been devised, primarily aimed at the measurement of dynamic interfacial tensions 



The pioneering theoretical work of Ward and Tordai [§] formulated a time- dependent 
relation between the surface density of surfactants adsorbed at an interface and their 
concentration at the sub-surface layer of solution, assuming a diffusive transport from the 
bulk solution. Consequent theoretical works have focused on providing a second closure 
relation between these two variables. Various relations have been suggested, resembling 
equilibrium isotherms 0, ffl, or having a kinetic differential form [5|, ffl. Such theories have 
been quite successful in describing the experimentally observea adsorption of common 
non-ionic surfactants. Yet, they suffer from several drawbacks: (i) The closure relation 
between the surface density and sub-surface concentration, which expresses the kinetics 
taking place just at the interface, is introduced as an external boundary condition, and 
does not uniquely arise from the model itself; (ii) The calculated dynamic surface tension 
relies on an equilibrium equation of state, and assumes that it also holds out of equilibrium 
0; (iii) Similar theories cannot be successfully extended to describe more complicated, 
ionic surfactant solutions ||. In the current paper we would like to briefly present an 
alternative approach, overcoming these three drawbacks. 

Consider an interface between an aqueous solution of non-ionic surfactants and an air 
or oil phase at x — 0. At x — > oo, the solution is in contact with a bulk reservoir of 
surfactants, where the chemical potential and surfactant volume fraction are fixed to be 
yu,;, and 0&, respectively. We consider a dilute solution, i.e. the surfactant volume fraction 
is much smaller than unity throughout the solution. The concentration is also smaller 
than the critical micelle concentration (cmc), so the surfactants are dissolved only as 
monomers. At the interface itself, however, the volume fraction may become large. 

We write the excess in free energy per unit area due to the interface (i.e. the change in 
interfacial tension) as a functional of the surfactant volume fraction in the bulk solution, 
4>(x > 0), and its value at the interface, (fi , Aj[(f>\ = J °° Af[(f>(x)]dx + /o(0o)- The first 
term is the contribution from the bulk solution, whereas the second is the contribution 
from the interface itself. The sharp, "step-like" profile considered has led us to treat 
the bulk solution and the interface as two coupled sub-systems, rather than a single one 
0. The bulk sub-system is considered as an ideal, dilute solution, including only ideal 
entropy of mixing and the contact with the reservoir, and neglecting gradient terms, 



where a denotes the surfactant molecular dimension and T the temperature (we set the 
Boltzmann constant to unity). However, since 0o may become much larger than <p(x > 0), 
we must take into account at the interface the finite molecular size and interactions 
between surfactant molecules, 



/o(0o) = TOo In 0o + (1 - 0o) ln(l - O )] - «0 O - (/V 2 )0o - ViM/a 2 - (2) 



The first term is due to the entropy of mixing, recalling that 0o is not necessarily small; the 
second accounts for the energetic preference of the surfactants to lie at the interface; the 
third is the energy of lateral attraction between neighboring surfactants at the interface; 
and the last term accounts for the contact with the solution adjacent to the interface, 
where the chemical potential is n\. 

Variation of A7 with respect to <fi(x) yields the excess in chemical potential at a 
distance x from the interface, A/i(x) = a 2 5A / y /5<f)(x) = T\n(j)(x) — fib for x > 0, and 
A/i = a 2 5Aj /5(f>o = Tln[0 o /(1 — <fio)) — a — j3(f>o — fi\. From these equations we can 
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deduce, as expected, //& = Tin 0& and fi% = Tln0i, where 0i denotes the sub-surface 
volume fraction. 

Thermodynamic Equilibrium. In equilibrium the chemical potential is equal to ^ 
throughout the entire system (the variations of A7 vanish). From the variation with 
respect to <j)(x) we obtain the equilibrium profile, <p(x) = 0& for x > 0, and from the 
variation with respect to 0o the equilibrium adsorption isotherm, 

0o = 6 /[0 6 + e- (Q+Wo)/T ]. (3) 

We have recovered the Frumkin adsorption isotherm, which reduces to the well-known 
Langmuir adsorption isotherm |10| when the interaction term is neglected ((3 = 0). From 
Eq. (§) and the above variations one also obtains the equilibrium equation of state, A7 = 
[T ln(l — O ) + (/?/2)0o]/ a2 ; expressing the equilibrium dependence of the surface pressure, 
II = — A7, on the surface coverage, 0o, according to a lattice-gas model. 
Out of Equilibrium. We assume proportionality between velocities and the chemical po- 



tential gradient [fLl"|, and take the surfactant mobility to be D/T, according to the Einstein 
relation (D being the surfactant diffusivity) . At positions not adjacent to the interface 
we obtain for the surfactant current density j(x) = —<f>(D/T)dfi/dx = —Ddcp/dx, and 
applying the continuity condition, d<p/dt = —dj/dx, get the ordinary diffusion equation, 
d<p/dt = Dd 2 <p/dx 2 . 

The proximity of the interface requires a more careful treatment. First, we discretize 
expression for A7 on a lattice with cells of length a, Aj[<fi) = a YhIi A/(0j) + fo((f>o), 
where <pi = <ft(x = ia). Discretized current densities, ji, can be similarly defined. Since 
we do not allow molecules to leave the interface towards the other phase (jo = 0), we 
have from continuity dfa/dt = —ji/a, and can therefore write dfajdt = — (j 2 — = 
(D/a)d(fi/dx\ x=a — d(po/dt. Applying the Laplace transform to the equations for d(f>/dt 
and d(f)\/dt while assuming an initial uniform state, <ft(x, t = 0) = <fib, a relation is obtained 
between the surface coverage and sub-surface volume fraction, 



Mt) = (l/a)y/D/^[2<j) b Vi- / h(r)(t - r)- l ' 2 dr] + 2<j) b - fa (4) 



This relation is similar to the classical result of Ward and Tordai @ , except for the term 
20b — 0i. The difference is due to fine details we have considered near the interface and 
our initial condition, and it vanishes when a goes to zero. Finally, we find the equation 
determining the kinetics at the interface itself, 

dfa/dt = faDQn - ^)/a 2 T = ( J D/a 2 )0 1 {ln[0 1 (l - O )/0 O ] + a/T + /?0 O /T} (5) 

Note, that the conditions at the interface are very different from those inside the aqueous 
solution, and the diffusivities, D, appearing in the equations above, cannot be expected 
to have strictly the same value. Solution of Eqs. (f|) and @ allows one to find the 
time-dependent surface coverage, <fio(t). 

By writing the above equations, we have separated the kinetics of the system into two 
coupled kinetic processes. The first takes place inside the bulk solution and is described by 
Eq. (|j), whereas the second takes place at the interface and is described by Eq. (|5|). Two 
limiting cases correspond to the relative speed of these two processes, (i) Diffusion-limited 
adsorption applies when the process inside the solution is much slower than the one at 
the interface. One can then assume that the interface is in constant equilibrium with 
the adjacent solution, i.e. the variation with respect to 0o vanishes, so 0o immediately 
responds to changes in X via the equilibrium isotherm, (ii) Kinetically limited adsorption 
takes place when the kinetic process at the interface is the slower one. In this case, the 
solution is assumed to be in constant equilibrium with the bulk reservoir, i. e. the variation 
with respect to <fi(x) vanishes [<p(x > 0) = 0J, and 0o changes with time according to 
Eq. (HI). 
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Looking at Eq. (^) we can identify the time scale of diffusion for attaining the equi- 
librium coverage, (po,eq, 

r d = (<po,e q /<p b ) 2 a 2 /D (6) 

Characteristic values of a 2 /D correspond to very short times (on the order of nanosec- 
onds), but the prefactor of (4>o, eq /4 ) b) 2 is typically very large (on the order of, say, 10 11 ). 
Thus, the diffusion time scales may reach minutes, as indeed observed in practice. In 
order to estimate the time scale of kinetic adsorption at the interface, we examine the 
asymptotic behavior of Eq (|5]) close to equilibrium and find [] 4>o teq — 4>o{t) ~ e~ t ^ Tk , with 

r k * (<W06) V/£)e- (a+Wo - )/T (7) 

Since the value of D at the interface is not expected to be drastically smaller than that 
inside the solution, comparison of Eqs. (H and (|7|) leads to the conclusion that > rfc. 
This result is somewhat expected, since we did not include any potential barrier in fo- 
Adding such a barrier, +e0o> to Eq. ([|) will result in a factor of e £ / T in 77, describing a 
kinetic limitation. We thus expect, in general, that non-ionic surfactants should exhibit 
diffusion-limited adsorption. This, indeed, has been observed for quite a large number 
of non- ionic surfactants JTJ. Q The "footprint" of diffusion- limited adsorption is the 

asymptotic time dependence |§] o , e g — 4>o(f) ~ \f r d/t- Any dependence between the 

surface tension and surface coverage will lead to a similar asymptotic time dependence of 
the dynamic surface tension as well. Four examples of experimental results are given in 
Fig. 1, all exhibiting the expected asymptotic t~ x l 2 behavior. 

We return now to the interfacial tension during the process of diffusion-limited adsorp- 
tion. As stated above, in this limit the interfacial contribution, /o(0o); is & t its minimum 
all the time, and A7 is given therefore by J °° Af[(j)(x)]dx + [T ln(l - O ) + (/3/2)^]/a 2 . If, 
in addition, we neglect the bulk contribution (recalling that it completely vanishes at equi- 
librium) we are left with the equilibrium equation of state. Therefore, relating the surface 
tension to surface coverage via the equilibrium equation of state approximately holds also 
out of equilibrium. Note, that this statement is valid only in the case of diffusion-limited 
adsorption, where fo is at its minimum during the whole process. The dependence of 
A7 on 0o, as defined by the equation of state, is shown in Fig. 2a. Since (fio constantly 
increases with time, we expect the time dependence of A7 to qualitatively resemble the 
curve depicted in Fig. 2b. This, indeed, resembles the curves found in experiments (e.g. 
0). Note the almost constant slope in the beginning of the process; the surface coverage 
significantly changes without a corresponding change in the surface tension. It is a result 
of the competition between the entropy and interaction terms appearing in the equation of 
state. The surface tension will start falling roughly when the second derivative of A7 with 
respect to <po changes sign (see Fig. 2a), i.e. when 1 — O ~ (P/T)^ 1 ^ 2 . As one examines 
surfactant solutions of increasing bulk concentrations, this surface coverage will be at- 
tained earlier along the process, and the initial plateau will shrink, until finally vanishing 
behind the finite experimental resolution. This trend is indeed observed experimentally 

E3- 

We have presented above an alternative model for the kinetics of non-ionic surfactant 
adsorption at fluid/fluid interfaces. It is a more complete model, in the sense that the 
kinetics in the entire system, in the bulk solution as well as at the interface, are described 
without the need for an additional, externally inserted boundary condition. We have 

1 Close to equilibrium we can also write Eq. (||) as d(f>o/dt~ (D(j)b/a 2 4>o,eq)[e ( - a+,3<l ' ^ T (f>i(l — </>o) — </>o], 
which coincides with the adsorption-desorption form of the Frumkin (or Langmuir, when /3 = 0) kinetic 
equation used by previous authors || ||. 

2 Note, that in the discussion above we have completely neglected a third time scale — that needed 
for lateral diffusion and molecular re-orientation at the interface. If, however, due to certain molecular 
constraints, this time scale is no longer negligible, exceptions to the above conclusions are to be expected 
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Figure 1: A variety of non- ionic surfactants were experimentally found to exhibit diffusion-limited 
adsorption. Four examples of dynamic surface tension measurements are given: 9.49 x 10 _5 M of decyl 
alcohol (open circles), from Ref. 0]; 2.32 x 1(T 5 M of Triton X-100 (squares), from Ref. f§; 6 x 1CT 5 M of 
Ci 2 E0 8 (triangles) and 4.35 x 1CT 4 M of Ci PY (solid circles), both from Ref. @. Note the asymptotic 
ir 1 / 2 behavior, characteristic of a diffusion-limited adsorption, and shown by the solid fitting lines. 



shown that relating the dynamic surface tension to surface coverage via the equilibrium 
equation of state, a procedure employed by practically all previous authors, is justified 
only in the case of diffusion-limited adsorption. Since the adsorption of non-ionic sur- 
factants is generally limited by diffusion, this assumption did not affect the validity of 
previous theories. We do not expect similar theories to be applicable to kinetically limited 
systems, such as salt-free ionic surfactant solutions ||. In contrast, our model allows, 
using the definition of A7 given above, for the calculation of dynamic surface tensions 
regardless of such limitations. Like any other free-energy formulation, the model can be 
easily extended to include additional interactions. A natural candidate should be the 
electrostatic interaction, i.e. extension of the model to ionic surfactants. This problem 
will be addressed in a future paper [pj|. 
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